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ABSTRACT 

Aims. As part of the WISH (Water In Star-forming regions with Herschel) key project, systematic observations of HtO transitions 
in young outflows are being carried out, with the aim of understanding the role of water in shock chemistry and its physical and 
dynamical properties. 

In this paper, we report on the observations of several ortho- and para-H20 lines performed with the HIFI instrument towards two 
bright shock spots (R4 and B2) along the outflow driven by the L1448 low-mass proto-stellar system, located in the Perseus cloud. 
These data are used to identify the physical conditions giving rise to the H2O emission and infer any dependence with velocity. 
Methods. We used a large velocity gradient (LVG) analysis to derive the main physical parameters of the emitting regions, namely 
n(H2), Tidn, A^(H20) and emitting-region size. Comparison has been made with other main shock tracers, such as CO, SiO and H2 and 
with shock models available in the literature. 

Results. These observations provide evidence that the observed water lines probe a warm (ryn ~ 400-600 K) and very dense (n ~10''- 
10' cm"^) gas, not traced by other molecules, such as low-y CO and SiO, but rather traced by mid-IR H2 emission. In particular, H2O 
shows strong differences with SiO in the excitation conditions and in the line profiles in the two observed shocked positions, pointing 
to chemical variations across the various velocity regimes and chemical evolution in the different shock spots. Physical and kinematical 
differences can be seen at the two shocked positions. At the R4 position, two velocity components with different excitation can be 
distinguished, with the component at higher velocity (R4-HV) being less extended and less dense than the low velocity component 
(R4-LV). H2O column densities of about 2 10'^ and 4 lO'"* cm"^ have been derived for the R4-LV and the R4-HV components, 
respectively. The conditions inferred for the B2 position are similar to those of the R4-HV component, with H2O column density in 
the range 10'"* - 5 10'"* cm"^, corresponding to H2O/H2 abundances in the range 0.5 - 1 10"'. The observed line ratios and the derived 
physical conditions seem to be more consistent with excitation in a low velocity I-type shock with large compression rather than in a 
stationary C-shock, although none of these stationary models seems able to reproduce all the characteristics of the observed emission. 

Key words. Stars: formation - Stars: low-mass - ISM: jets and outflows - ISM: individual objects: L1448 - ISM: molecules 

1 . Introduction since its abundance can vary by order of magnitudes through the 

shock lifetime (e.g. Bergin et al. 1998). Due to the difficulty of 

Strong radiative interstellar shocks are produced by the inter- observing water with ground-based facilities, the study of wa- 

action of supersonic mass ejections from young stellar objects ter is restricted to IR/sub-mm telescopes from space. Dedicated 

H with the dense ambient cloud. The observational signature of sub-mm satelfites, such as SWAS and Odin, have allowed, for 

^ these shocks are the bright line emissions from molecules and t^e first time, to observe the ground state ortho-HzO at 557 GHz 

atoms, whose excitation conditions and relative abundances re- and t o compare its profile and abundanc e with that of CO (e.g. 

veal fundamental information on the type of interaction, and on Fra nkUn et al.ll2008[ lBierkeliltIDl20U9h . The poor spatial res- 

the physical properties of both the jets and the ambient medium, ojution of these facifities, together with the restriction of ob- 

In the dense environments of young proto-stars (the so-called serving a single line, have however prevented the study of the 

Class sources) the gas cooling occurs through emission of Hz, water excitation conditions and their variations along outflows. 

CO and H2O over a wide range of wavelengths, spanning from Through the Infrared Space Observatory (ISO), on the other 

near-IR to sub-mm wavelengths (e.g. Kaufman & Neufeld 1996; hand, it has been possible to infer the physical conditions of 

[Flower & Pineau Des Forets 2010). Among these main coolants the warm co mponent giving rise to stro ng water emission (e.g. 

of outflow shocks, water is the most sensitive to local conditions, iLiseau et al.l [i996: Giannini etal 1 12001 . Nonetheless, the lim- 

ited spatial and spectral resolution have not allowed the associ- 

* Herschel is an ESA space observatory with science instruments ation of this warm gas with a specific kinematical component, 

provided by European-led Principal Investigator consortia and with im- The Herschel Space Observatory is now able to overcome these 
portant participation from NASA. 
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limitations, thanks to its improved spatial and spectral resolution 
over a large wavelength range. 

As pa rt of the WISH (Water I n Star-forming regions with 
Herschel, Ivan Dishoeck et al ] l20Tlh kev program, we are under- 
taking systematic observations of H2O transitions in young out- 
flows, employing the two spectrometers on board Herschel, HIFI 
and PACS. One of our main targets of investigation is the L1448- 
mm outflow. This flow, powered by a low-luminosity ( 11 Lp) 
Class proto-star located in Perseus {d - 235 pc, Hir ota et alJ 
Hon]), is known to be a strong water emi tter on the basis of 
ISO observations (Nisi ni et al.iri999l l200dh . HIFI observations 
of water in the central region of the LI 448 outflow (within 20" 
radius from the source) have been presented bv Kristensen et aTl 
(l20Tll) . who focused on the H2O properties in the extreme high 
velocity gas associated with the source collimated jet. 

In this paper, we present HIFI observations of several H2O 
transitions obtained in two active shocked regions along the 
L1448 flow. Our aim is to study the excitation conditions of H2O 
in the shocked gas, exploring, in particular, variations of excita- 
tion with velocity. We want in this way to study the physical and 
chemical conditions of the interstellar medium, which has been 
affected by the shock. 



2. Observations and data reduction 

Figure [U shows the JCMT CO(3-2) emission in contours, over- 
laid on the PACS map of the H2O 2i2-loi emission at 1670 GHz 
towards L1448 (Nisini et al. 2011, in prep.). The H2O map ex- 
hibits several emission peaks, which correspond to positions of 
active shock regions, named B2-B3 (for the blue-shifted lobe) 
and R2-R3-R4 (for the red-shifted lobe), following the nomen- 
clature of Bachilleret al. (1990), based on the analysis of CO 
line profiles. For our water line survey we selected the shock 
spots B2 an d R4: the first corre sponds to a bright H2 emis- 
sion region (lDavis&Smiflilll995h . arising from the interaction 
of the high velocity jet from the mm source with the ambi- 
ent medium; while the second is located at the end of the red- 
shifted lobe and is spatia l ly asso ciated with the bow-shock seen 
in SiO by iDutrev et al.l (Il997l) . Their offsets with respect to 
flie central driving source L1448-mm (q'(J2000)=03''25™38!9, 
(J(J2000)=-h30°44'05'.'4) are respectively (-16", 34") forB2 and 
(26", -128") for R4. 

A survey of several ortho- and para-H20 li nes has been con- 
ducte d with the HIFI heterodyne instrument (de Graauw et alJ 
2010| ) on board the Herschel Space Observatory (Pilbrattet al. 
2010 ). towards these two positions. The survey comprises lines 
with excitation energies ranging from 27 to 215 K and also in- 



30°46'00" - 



eludes the ortho-Hj^O lio - loi transition at 547.7 GHz, useful 
to infer opacity effects. 

The data were processed with the ESA-supported package 
HIPlJ] {Herschel Interactive Processing Environment, Ott 2010) 
for calibration. The calibration uncertainty is taken to be 20%. 
Severe baseline problems have been found at the B2 position 
in Band 6 and Band 4, corresponding to the ortho-H20 2i2 - loi 
line at 1670 GHz and the 3i2-3o3 line at 1097 GHz, respectively. 
We obtained, using the experimental "matching technique'^ for 
electronic standing waves, a tentative detection for the former 



' HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 

^ The techniques are being developed at the HIFI ICC; see 
http://herschel.esac. esa.int/twiki/bin/view/ I^blic/HifiCalibrationl 
Web?template=viewprint for references. 
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Fig. 1: PACS image of L1448 at 1670 GHz in false colors, with 
the negative contours in dotted line, and the JCMT CO(3-2) 
emission (Nisini et al. 2011, in prep.) in contours: the blue- 
shifted emission is integrated between -100 and 4 km s"' and 
the red-shifted emission between 6 and 100 km s The central 
position of the map, as well as the positions chosen for the HIFI 
line survey (R4 and B2), are marked with a green cross. For the 
latter two positions the largest and the smallest HIFI beam sizes 
are shown in green. 



line at a ~3cr level, which is consistent with the PACS map at 
1670 GHz, presented in Fig.[TJ for the latter instead only an up- 
per limit has been derived. 

Further reduction of all the spectra, including baseline sub- 
traction, and the analysis of the data were performed using the 
GILDA^ software. H- and V-polarizations were co-added af- 
ter inspection; no significant differences were found between 
the two data sets. The calibrated T* scale from the telescope 
was converted into the r,nb scale using the main-beam efficiency 
factors provided by Roelfsema et al. (201 1, submitted^ and re- 
ported in Table[T] The beam sizes range from ~13" to 39". The 
largest and the smallest beam sizes of the HIFI observations are 
marked as green circles in Fig. [T] for each observed position. 
At the velocity resolution of 1 km s"', the rms noise range is 
3-100 mK(r,nb scale). 

Transitions from other molecules have also been observed in 
the HIFI bands centered on the water lines, namely NH3(lo-Oo), 
i^CO(10 - 9), HCO+(6 - 5) and C"*0(5 - 4). The only detected 
line among these is the NH3 transition at the B2 position, while 



3 http://www.iram.fr/IRAMFR/GILDAS/ 

* see also http://herschel .esac.esa.int/twiki/bin/view/Public/| 



HifiCalibrationWeb?template=viewprint 
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Table 1 : Summary of the HIFI line survey. 
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" Integrated intensities of the lines in the two velocity components at R4 (between and 20 km s"' for R4-LV and between 20 and 60 km s"' for 
R4-HV) and total integrated intensities of the lines at B2. 
Total observing time: on + off + overheads. 
All spectra were smoothed to a velocity resolution of 1 km s"'. 

Upper limits to the integrated intensities, computed from the 3cr values of the spectra and the line widths of the respective velocity components. 
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Fig. 2: HIFI spectra at the R4 (left) and B2 (right) positions. Note the different profiles for low-lying (black) and high-lying (red) 
levels at the R4 position. 



for all the other transitions we have non-detections at a rms noise 
level (at the velocity resolution of 1 km s ') of 4 mK (HCO^ 6-5 
and C"^0 5-4) and 20 - 26 mK ('^^CO 10 - 9). A summary of 
the observations, including frequencies (v), beam sizes (HPBW), 
main-beam efficiencies [7]^^) and velocity-integrated intensities 
is given in Table [T] 



3. Line profiles 

The spectra of the H2O lines observed at the R4 and B2 spots are 
presented in Fig.|2] respectively on the left and right panels. 



At the R4 position, the water line profiles reveal a broad 
emission extending up to about 50 km s"' with respect to 
the systemic velocity. The overlay between water line profiles 
with different upper level energies shows that the lines having 
iiu < 137 K peak at velocities around 25 km s"', while the lines 
with higher upper level energies peak at lower velocities (around 
10 km s '). The diff'erent profiles of lines at diff'erent excitation 
could in principle be caused by large self-absorptions in the low 
velocity range of those lines connecting with the ground state. 
This could be due to the presence of large columns of cold wa- 
ter along the line of sight. However, we do not think this is the 
case, since the effect is not observed in the spectra obtained at 
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Fig. 3: Overlay between the HIFI H2C ) 557 GHZ li ne (bl ack), 
IRAM-30m SiO(2-l) line (red) from iNisini etall (l2007h and 
JCMT CO(3-2) line (green) from Nisini et al. (2011, in prep.), 
towards the R4 (upper) and the B2 (lower) positions. 

the B2 position. In addition, as we will show later in Fig. [3] 
the profile of the H2O lio - loi line at 557 GHz matches well 
the SiO 7=2-1 line profile, which does not show any evidence 
of self-absorption even at the systemic velocity. The two veloc- 
ity regimes at R4 (and correspondingly the different excitation 
lines) are thus probably tracing gas with different physical con- 
ditions. 

In particular, from the comparison of the shape of the profiles 
shown in Fig.|2l it can be noted that the lines at higher excitation 
(namely the 0-H2O at 1097 GHz and the P-H2O at 752 GHz) 
have a simple triangular profile peaking at low velocity, simi- 
lar to those observed also in B2, while the other lines present 
an additional component peaking at higher velocity, overlaid on 
the low velocity triangular profile. It thus seems that the ob- 
served profiles results from the superposition of two gas com- 
ponents, having different kinematical and excitation properties, 
seen along the line of sight. 

Finally, the non-detection of the ortho-Hj^O lio - loi line 
at 548 GHz is shown in the bottom-left panel of Fig. |2l with a 
rms noise of about 3 mK (T^b scale) at 1 km s"' velocity res- 
olution. This gives a constraint on the H2O lio - loi optical 
depth. Assuming an oxygen isotope ratio '^0/"^0~ 560 from 
[Wilson & Rood ( 1994, for local interstellar medium) and under 
the assumption that the two lines have the same Tex, we find that 
T557GHZ < 15 at the peak of the H2O emission (based on Hj^O 
3cr upper limit). 

In contrast with R4, all the detected lines at the B2 position 
show a similar triangular line profile with a line wing extending 
up to about -50 km s ' . It is interesting to point out that, among 
the detected H2O lines, the highest excitation line (2ii-2o2 with 
iiu = 137 K) seems to show emission extending at high velocity, 
around -45 km s"\ at a level of < 3 cr. This may indicate the 
presence of high excitation gas at high velocity. 

The bottom-right panel of Fig.|2]also shows the NH3(lo-Oo) 
emission line detected at B2. Ammonia emission is confined at 
the systemic velocity of about 5 km s"', testifying to its ori- 
gin from the molecu lar cloud and not from the outflow (see also 
liBachineretaDl990h . 



Figure [3] shows the comparison between the H2O lio - loi 
line profile in black and that of SiO v=0, 7=2-1 in red and the 
CO(3-2) in green (Nisini et al. 2007, 201 1), towards both investi- 
gated positions. The figure highlights the very different line pro- 
file behaviour at the two positions. In particular, the profile of 
the H2O line at 557 GHz follows quite closely the SiO profile 
at the R4 position, while at the B2 position the SiO emission is 
enhanced at the extreme high velocities (EHVs, v ~ 55 km s"'). 
Conversely, there is no water enhancement above the wing pro- 
file at the highest velocities, which suggests that the water emis- 
sion from the EHV gas is observed only closer to the driv- 
ing source (see also K ristensen et ani201 I h Nisini et al. 2011, 
in prep.). A discussion of the different behaviour of H2O and 
SiO line profiles at the two observed positions is given later, in 
Sect.|5] 

At both positions, the comparison with the CO(3-2) emission 
points out that H2O always has excess emission with respect to 
CO at intermediate velocities, suggesting that H2O and CO could 
trace different gas components. This will be further discussed in 
Sect. 121 

4. Excitation analysis 

The physical conditions of the observed emitting regions have 
been investigated by comparing the observed intensities with 
predictions from t he RADEX esc ape probability code in plane 
parallel geometry jvan der Tak et al. 2007). Although the plane 
parallel geometry i s not able to reproduce the line profiles (e.g. 
Bj erkeli et alj |201 1), it can be used as a first approximation 
for this study. The molecular data were taken from the Leiden 
Atomic and Molecular Database (LAMBDAB Schoier et alJ 
2005), with the collisional rate coefficients from Faure et aL 
(120071) . 

As pointed out in Sect. [3] the profiles at the R4 position 
present two distinct components showing a clear change in ex- 
citation with velocity. We therefore need to separate the two ve- 
locity components and analyze them separately. An approach to 
separate the two velocity components in each line may be to fit 
the triangular component, peaking at low velocity, in each water 
line and subtract it from the rest of the line. However, due to the 
fact that fitting a triangular shape to the water lines is arbitrary, 
we opted for a simpler approach. We divided the emission at R4 
into two velocity components: a low-velocity component (R4- 
LV hereafter, between and 20 km s"') and a high-velocity one 
(R4-HV hereafter, between 20 and 60 km s~^). As described in 
the Appendix, the results of the analysis will be little affected 
by the method we used to separate the two velocity contribu- 
tions, due to the fact that both components dominate the relative 
velocity range we considered. 

For B2, instead, due to the low signal-to-noise and the non- 
detection of one H2O line out of five, we had too little informa- 
tion to perform an excitation analysis as a function of velocity, 
and therefore we have integrated the emission over the whole 
velocity range. The observed velocity-integrated intensities for 
each component are given in Table[T] 

We have compared the integrated intensities with a grid of 
RADEX models, constructed by varying the parameters in the 
following ranges: n(H2) = lO'* - 10** cm ^ A?(ortho-H20)= 
IO12 - lO'*" cm-2, rkin=100-1600 K and ortho to para ratio 
o//? = 3, 2.5, 2, 1.5. Line widths (FWZI) of 20, 40 and 50 km s"' 
have been adopted for R4-LV, R4-HV and B2, respectively. 
Given the different beam sizes of the observations, we have also 

' http://www.strw.leidenuniv.nl/ moldata/ 
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considered the size of the emission region (0) as an additional 
parameter: we varied it from almost point-like (few arcsec) to 
extended emission. A Gaussian emitting region, as well as a 
Gaussian beam, were assumed to correct the emitting size for 
beam dilution effects. 

4.1. R4 position 

The comparison of the observations with the RADEX grid of 
models shows that only a limited range of physical parameters 
is able to simultaneously reproduce both the line ratios and the 
line intensities in each component (see Figs. IA.l] and lA.2l in the 
Appendix). In Fig. |4] these best fit models are compared with 
the measured integrated intensities. For the R4-LV component, 
we are able to reproduce the observations, within the calibra- 
tion errors, assuming only extended warm gas (Tyn ~ 600 K) 
having a very high density n(H2) ~ 10^ cm"^ (model LV- 
1 in Table |2]i. The implied column density is of the order of 
A^(H20) ~ 2 10^^ cm"^, which indicates that the observed lines 
have low opacities (tssvghz < 0.1), consistent with the non- 
detection of the Hj^O transition. 

For the R4-HV component, on the other hand, there is more 
degeneracy in the physical parameters (see Fig. IA.2I ) and the 
observations can be reproduced either by 10^ cm"-' gas at low 
temperature (Tyn ~ 150 K, model HV-1) or by gas as warm 
as in the R4-LV component, but having a lower density of 
~ 10* cm"^ (model HV-2). It is therefore not possible, at this 
stage, to infer whether the lower excitation of the R4-HV com- 
ponent with respect to the R4-LV component, evidenced by the 
comparison of the line profiles, is due to a temperature or a den- 
sity effect. In both cases, the R4-HV component appears to be 
more compact (~ 13-21") and with a higher column density 
A^(H20) = (7 - 40) 10'^^ cm^^ with respect to the R4-LV gas. 

The different sizes found in the R4-LV and R4-HV compo- 
nents are consistent with the PACS map at 1670 GHz presented 
in Fig.[T] The map in fact shows that the 1670 GHz emission at 
R4 is dominated by a compact emission, which may correspond 
to R4-HV, superimposed on a more extended weak emission, 
which may correspond to R4-LV. 

The physical conditions derived for the two components, in 
particular the very high density, can be compared with those in- 
ferred from other shock tracers, to see whether the observed H2O 
emission probes the same gas. Given the observed similarity be- 
tween the H2O 557 GHz and the SiO(2-l) line profiles shown in 
Fig. |3] one would be tempted to conclude that the two species 
have similar excitation conditions. To test this hypothesis, we 
have used the SiO multi-line observations at R4 presented by 
iNisini et~al ](|2007) to derive temperature and H2 density of the 
SiO gas in both the R4-LV and R4-HV components. We have 
used a grid of LVG models, constructed from the RADEX code, 
and assumed an emission size of the order of 15", for both com- 
ponents, estimated from the IRAM Plateau de Bure (PdB) chan- 
nel m aps of the SiO(2-l) emission, presented by iDutrev et al.l 
(fT997.) . 

The best fit models that reproduce the SiO emission in the 
two velocity components are given in Table |2](LVsio-2 for R4- 
LV and HVsio-3 and HVsio-4 for R4-HV). The models trace 
warm gas (T\^n between 250 and 600 K) but with H2 densities 
of only (1 - 5) 10"^ cm"^, thus two order of magnitudes lower 
than those derived from the H2O analysis. In particular , the H 2 
densities are lower than those derived from'Ni sini et a n (I2007h . 
n(H2) ~ 2.5 10^ cm"^ with T = 200 K, and this is probably 
due to beam-filhng effects that were not taken into account in 
iNisini et"an (l2007h They assumed in fact a beam-filling factor 



Table 2: Summary of the models derived for each component. 
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1 D 


-^T- 




A'(H20)'' 











y^^j 






(tircscc) 


R4-LV 


LV-l" 


3 


600 


10' 


210" 


37 




LVsio-2''' 


3 


500 


3 10-* 


3 lO'"' 


15 


R4-HV 


HV-1 


3 


150 


3 10' 


7 10" 


21 




HV-2'' 


3 


650 


10^ 


410'" 


13 




HVsio-3'^ 


3 


250 


5 10* 


9 10'5 


13 




HVsio-4'^^ 


3 


600 


10^ 


3 10"^ 


13 


B2 


B2-1 


3 


450 


6 10^ 


10'* 


33 




B2-2''-'' 


3 


450 


lO*^ 


5 10'* 


17 



Ortho- + para-H20. 
'' The bold font highlights for each component the model that better 

fits all the observational constraints. 
^ Physical parameters (r^n, nH,'®) derived from the SiO analysis 

(see Sect.l4ll. 
'' See Sect.l42l 



7 R4-LV 




r 


1 ' 

I 


1 1 1 1 _ 

1 -7 

Model LV-1 - 
Model LVsio-2 _ 


= — — ^ — ^ 


T 












i 








Model HV-1 - 










Model HV-2 - 










Model HVsio-3 - 


7 R4-HV 








Model HVsio-4 _ 


= — — ^ — ^ 
7 B2 


1 1 

I 




f 

1 , 


Model B2-1 - 
Model B2-2 _ 

1 1 1 1 - 



500 1000 1500 

Frequency (GHz) 



Fig. 4: Top and middle panels: Comparison between the ob- 
served water intensities and the different models in Table |2] for 
R4-LV (top), R4-HV (middle) and B2 (bottom). The intensities 
predicted by the models have been coiTected for the relative pre- 
dicted filling factors. The given errors include 20% of calibration 
error. The triangles represent the models derived from H2O and 
the squares from the SiO. Bottom panel: Same comparison as the 
top panels, but for B2. The upper limit of the non-detected H2O 
line at 1097 GHz is also shown for each model. 



equal to 1, which implies that the line intensities are reproduced 
with a lower column density and consequently a higher particle 
density. 

We have then investigated how well the conditions estimated 
from the SiO emission reproduce the water lines, by varying 
the H2O column density in order to match the 557 GHz H2O 
line. The results are visualized in Fig. |4] For R4-LV, the SiO 
best fit model underestimates the intensities of the water lines at 
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higher excitation (para-HaO 2ii -2o2 at 752 GHz and ortho-H20 
3i2 -3o3 at 1097 GHz). This could indicate that SiO is tracing an 
additional low-velocity gas component at lower excitation than 
that probed by H2O. Finally, for R4-HV we considered two pos- 
sible SiO models (HVsio-3 and HVsio-4 in Table |2]i and both of 
them reproduce reasonably well only four out of the five H2O 
lines. 

The high density regime found by the best-fit models, based 
only on the water lines, are more consistent with the conditions 
deriv ed from Spitzer mid- IR H2 observations along the L1448 
flow (iGiannini et al.ll20TTI) . Although the R4 position is not cov- 
ered by these observations, the other H2 shocked spots in the 
red-shifted lobe are consistent with a warm gas at densities of 
the order of 10^ cm"^. 

4.2. B2 position 

At the B2 position we detect four (the 1670 GHz water line is 
however only a tentative detection) out of the five targeted tran- 
sitions: therefore we have too little information to constrain all 
the parameters of the fit (see Fig. lA.3l in the Appendix). In order 
to overcome this problem, we assume that the water emission 
originates from the same gas emitting the mid-IR H2 lines, ob- 
served with Spitzer (Giannini et al. 2011). This assumption is 
supported by the results found in other out flows observed within 
the WISH program (e.g. lNisini et alJ2010l Tafafla et al. 201 1, in 
prep.). The H2 rotational transitions at low-7 indicate a tempera- 
ture of about 450 K, while a non-LTE fit through all the H2 lines 
suggests high densities, in excess of lO*" cm"-'. 

Keeping these constraints, the fit is still degenerate in the 
product n(H2)xA^(H20), i.e. the data can be reproduced either 
by n(H2)~6 lO*- cm^^ and A^(H20)= 10"* cm^^ (model B2- 
1 in Table |2] and Fig. |4| or n(H2)~10^ cm ^ and A^(H20)= 
5 10'** cm"^ (model B2-2). The smaller size of the B2-2 model, 
with respect to the B2-1 model, seems to be consistent with the 
rather compact PACS emission at the B2 position, as shown in 
Fig. [U Hence, the B2-2 model is indicated in Table |2] as the 
best model for the B2 component, even if it is associated with 
a worse fit. In fact, the comparison between the integrated inten- 
sities, predicted by the B2-2 model, as seen in FiglH shows that 
the intensities of the H2O lines with the biggest beam sizes (at 
557 GHz and 752 GHz) are slightly underestimated. This could 
be due to the fact that the observations performed with these 
large beam sizes may collect emission from additional compo- 
nents, especially from the source position, as can be evidenced 
from the PACS map of Fig. □ 



5. Discussion 

5.1. Water abundance 

The presented excitation analysis shows that the observed HIFI 
water lines are consistent with warm gas at high density hav- 
ing moderate H2O column densities that do not exceed ~5 10'"* 
cm"^. Similar water conditions have also been found in the 
LI 157 flow (Vastaetal. 2011). These water column densities 
are significant ly lower than those der ived close to the L1448- 
mm source by iKristensen et al.l (|20I li) . who estimated A^(H20) 
^ 10'^ cm"^ in all the considered velocity co mponents. This is 
in agreement with ISO observations of L1448 (Nisini et al .l 19991 
[2OOO), where the strongest water emission is found towards the 
central source. Moreover, this seems to suggest that the water 
emission from the EHV gas is only observed at high column den- 



sities, close to the driving source (see IKristensen et ani201 ih . in 
contrast with other shock tracers. 

H2O/CO abundance ratios found bv IKristensen et all (1201 ll) 
are in the range 0.01-1, under the assumption that the two 
molecules are probing the same gas component. Adopting the 
same assumption, we have applied the physical conditions de- 
rived from water to the CO(3-2) data and computed the cor- 
responding CO column densities using the RADEX code. The 
data were corrected for beam dilution effects and the same sizes 
as estimated from water were assumed. The derived N(CO) val- 
ues are high in comparison with water, ranging between ~ 5 10'* 
at R4-LV and HV, and 2 10'^ cm ^ at B2. These values would 
imply a H2O/CO abundance ~ 10 "'-lO"^, if one assumes the typ- 
ical CO/H2 ratio of 10 i.e. a water abundance extremely low 
for a shocked region. However the assumption that the H2O and 
CO(3-2) emission come from the same region is likely not cor- 
rect, given that they largely differ in both critical density (a factor 
of 10^) and excitation temperature, and given also their different 
line profiles (Fig. |3]l. It is thus conceivable that CO(3-2) traces 
mostly entrained ambient gas, while the H2O emission comes 
directly from the warm shocked gas. 

A more direct estimate of the H2O abundance of the gas com- 
ponent traced by the HIFI lines can be obtained comparing the 
water column density with the H2 column densi ty of the same 
component, inferred from Spitzer observations dGiannini et al.l 
2011; Dionatos et al. 2009). Only the B2 position was covered 
by these observations: the H2 gas is stratified in temperature, 
with the gas at T~ 400-500 K having A^(H2) of the order of 
5 10'^ cm"^ in a ~20" beam. The H2O abundance at B2, rel- 
ative to this gas component, can therefore be constrained to be 
~ 0.5 - 1 10"^. A similar abundance can be inferred for the 
R4-HV gas (whose estimat ed density is consist ent with the av- 
erage H2 density given bv iGiannini et alJl201 ll along the flow), 
assuming that the same A^(H2) applies to this component. This is 
clearly just a rough estimate, since we are comparing gas com- 
ponents at slightly different temperatures. These abundances are 
lower than those expected in non-dissociative, stationary shocks, 
where most of the gas-phase oxygen is converted into water 
( Kaufman & Neufeld 1996). The derived H2O abundances are 
however higher than interstellar values, indicating that the ob- 
served warm gas has been processed by shock-related chemical 
processes. 

5.2. Comparison witli shock modeis 

In order to better investigate the consistency of the observed 
emission with shock models, we have compared the observed 
intensity ratio s with the grid of stationary C- an d J-type shocks 
provided by Flower & Pineau Pes ForetsI (1201 01) . This grid ex- 
plores diff'erent shock velocities (from 10 to 40 km s~*) and two 
pre-shock densities (2 lO"* and 2 10^ cm^^). For each type of 
shock (C- and J-), we have identified the model that better re- 
produces our observed line ratios. As in the case of the RADEX 
analysis, the emission size has been taken as an additional pa- 
rameter, used to correct the intensity of lines with different beam 
sizes. The results are shown in Fig. |5] where the observed and 
predicted line ratios are plotted with respect to the 1097 GHz line 
for R4-LV and R4-HV and to the 1 1 13 GHz Une for B2. From 
the figure one can see that the line ratios at R4-HV are better 
reproduced by a J-type shock with a shock speed of 30 km s 
consistent with the line width of this velocity component. The 
shock conditions of the R4-LV component are more degenerate: 
in fact, both shock models seem to reproduce the data within the 
uncertainties. There are however different reasons to prefer the 
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Fig. 5: Comparison between the observed line ratios (circles) of the H2O lines and the corresponding line ratios predicted (triangles) 
by the best-fit C-type shock (top) and the J-type shock (bottom) models. For each type of shock (C- and J-), the model that best 
reproduces our observed line ratios is presented and the corresponding shock velocity, pre-shock density and emission size are 
marked. The line ratios are normalized to the 1097 GHz line for R4-LV (left) and R4-HV (middle) and to the 1113 GHz line for 
B2 (right). For B2, the detected lines are shown as black filled circles, while the black open circle represents the upper limit of the 
non-detected H2O line at 1097 GHz (assuming a 3cr upper limit and a line width of 50 km s '). 



J-type shock model. First of all, the J-shock model has a shock 
speed of 20 km s"' (versus 40 km s ' for the C-shock), consis- 
tent with the line width of the R4-LV component. In addition, the 
best-fit J-shock model have a pre-shock density of 2 10^ cm"^, 
higher than that predicted for R4-HV (2 10"^ cm ■'), and this is 
consistent with our LVG analysis, showing that the difference in 
excitation between the R4-LV and R4-HV components can be 
due to a different pre-shock density. The comparison between 
the post-shock density derived from the LVG analysis and the 
pre-shock density predicted by the J-shock, indicates a compres- 
sion factor npost/npie ~100, which is consistent with water emis- 
sion originating from the post-shock gas regions further behind 
the shock front, where H2, which has been initially dissociated, 
starts to reform, ke eping the temperature to values a round few 
hundreds of Kelvin jFlower & Pineau Pes Foretsll2010f) . 

At the B2 position, only the J-type shock predicts a shock 
velocity consistent with the observed wide line profiles (i.e. 
30 km s ' against 10 km s"' for the C-type shock), although both 
C- and J-type shock models reproduce the line ratios. This best- 
fit J-type model predicts a shock speed of 30 km s ' and a pre- 
shock density of 2 10^ cm"^. Several additional lines of evidence 
indicate that a dissociative shock must be present at this posi- 
tion. Spitzer-IRS observations show that atomic emission from 
Fell] 26 jum and [S i] 25 nm is associat ed with this position 
Neufeldet all 120091: bionatos et a"ni2009l) . R\CS observations 
reveal emission from OH and [Oi] 63//m (Santangelo et al., in 
prep.): this is a further indication that the gas-phase reactions 
that convert all the oxygen into water, activated at T > 300 K, 
have not been completed a nd dissociative shocks must be present 
to have high OH emission (iNeufeld & Dalgarnolll989h . 

There are however also several inconsistencies between the 
considered models and our observations: in particular, none of 



these models (either C- or J-type) is able to reproduce the abso- 
lute intensities of the observed lines. Indeed, the comparison be- 
tween the predicted and observed intensities would imply emis- 
sion sizes much smaller (of the order of few arcsec) than those 
fitted through the line ratios. We believe that these discrepancies 
are due to the simplified nature of the adopted shock models. 
In fact, a single, plane-parallel shock model is unlikely to pro- 
vide a completely satisfactory representation of the line emis- 
sion from a complex shock structure, where both geometrical 
and temporal effects can significantly alter the emergent line 
intensity and the final chemical abundance. In particular, sev- 
eral studies have shown that shocks from young outflows have 
not reached the steady state yet, and thus are better represented 
with non-stationary shock co nditions that combine together both 
C- and J-type shock waves jFlower et al.ll2003t iGiannini et al.l 
2006; Gusdorfetal. 2008b). Predictions of H2O hne intensi- 
ties in CJ-shocks are provided for a limited set of conditions by 
iGusdorf et al.l (1201 Ih : a qualitative comparison indicates that ab- 
solute intensities and water abundances derived in these models 
are similar to those observed in L1448. Clearly, the development 
of more specific models, suited for the shock regions under in- 
vestigation, would be needed for a better understanding of the 
water properties highlighted by our observations. 

5.3. Comparison with SiO emission 

Finally, we discuss the very different line profile behaviour of 
H2O and SiO at the two investigated positions (see Fig. [3] in 
Sect. ID. At R4, the profile of the H2O line at 557 GHz follows 
quite closely the SiO profile. The SiO abundance is enhanced 
in shocks due to both sputtering of silicon from the grain man- 
tles, efficient at shock velocities as low as ~ 10 km s"\ and core 
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grain disruption occurr ing at higher shock speeds, of the order 
of 25 - 30 km s ' (e.g. Jimenez -Serra et al1l2008l : lOusdorf et al.l 
l2008ab . The observed SiO emission in both the R4-LV and R4- 
HV components is therefore consistent with the range of shock 
speeds we infer from the H2O analysis. 

At B2, the comparison between water and SiO shows a more 
complex situation, with a strong variation of the SiO/H20 ra- 
tio as a function of velocity. This behaviour probably reflects 
very diff'erent SiO/H20 abundance ratios in the various velocity 
regimes, indicative of a different chemistry rather than a differ- 
ence in the excitation conditions. In particular, the SiO emission 
at B2 is enhanced at the EHVs (v ~ 55 km s"'), where no excess 
of water emission is detected. This is in contrast with the ob- 
servations of water towards the dri ving source, where the EHVs 
water emission is clearly detected (iKristensen et al.ll201 ll; Nisini 
et al. 2011, in prep.), which indicates that EHV water emis- 
sion is more localized than SiO and show up only in the regions 
with high water column densities. It has been suggested that the 
EHV gas traces the primary proto-stellar jet, whe re molecules 
are rapidly synthesized from the initial atomic gas ( Tafalla et al.l 
[2OIO). Available chemical models for the chemistry in the pri- 
mary jet, although simplified in the geometry and physical struc- 
ture, predict a low H2 0/SiO ratio for values of mass loss rates 
less than 10"^ M@ yr"' jGlassgold et al.ll99Ih . which is in agree- 
ment with the molecular mass flux rate of 10"^ yr"', derived 
forL1448-mmby Dionatos et al. (2009). Regarding the very low 
SiO associated with the wing emission, this is a bit puzzling 
given that these wings extend up to 50 km s i.e. where one 
would expec t sputtering to be ve ry efficient in releasing silicon 
from grains. Tafalla et al.l (|2010'), studying the chemistry of the 
L1448-R2 shock spot, located in the red-shifted lobe symmet- 
rically to the B2 position, point out a strong trend of increasing 
SiO abundance with velocity, going from the lower velocities up 
to the EHV gas, which is not followed by the other molecules 
they observed. It therefore appears that there are strong chem- 
istry variations across the various velocity regimes at different 
positions that are not easily explained by current shock chemical 
models. 

6. Conclusions 

Herschel observations of several ortho- and para-H20 transitions 
towards the outflow driven by the L1448 low-mass proto-stellar 
system have been presented. These observations are part of the 
WISH key program. Two shocked positions, R4 and B2, along 
the main L1448 outflow have been studied. The main results of 
this work are the following: 

1 . The two investigated positions (R4 and B2) show physical 
and chemical differences. In particular, R4 exhibits strong 
variations in the excitation conditions as a function of the 
LSR velocity and two velocity components can be distin- 
guished in the water emission: a low-velocity component 
(R4-LV) and a high-velocity component (R4-HV). The exci- 
tation at R4 decreases with velocity, i.e. the lines with higher 
upper level energies peak at lower LSR velocity, while the 
lines with lower upper level energies peak at higher LSR ve- 
locity. 

2. The observed emission in both shocked positions is best rep- 
resented by very dense (hh, ~ 10^-10^ cm"-') and warm 
(Tkin =400-600 K) gas, having moderate H2O column densi- 
ties: NiHiO)- 10'^ cm-^ at R4-LV and ~ 10'^-5 lO'^ cm"- 
at R4-HV and B2. These column density values correspond 
to H2O/H2 -0.5-1 10"^ at B2 and R4-HV, obtained from 



comparison with A^(H 2) derived from Spitzer observations 
(iGiannini et al.ll201ll) . 

3. The inferred physical conditions seem to be better repro- 
duced by a J-type shock, where large compression factors 
leading to high density in the post-shock gas are expected. 
Moreover, the relatively low observed column densities and 
the PACS OH and [OI] 63/im observations are consistent 
with models where H2O takes time to reform from fully 
dissociated gas. However, more detailed shock models are 
needed. 

4. These Herschel H2O observations provide evidence that wa- 
ter is unique in tracing very dense shock components, not 
traced by other molecules, such as CO and SiO, but rather 
traced by H2. In particular, H2O shows strong differences 
with SiO in the excitation conditions and in the line profiles 
in the two observed shocked positions, pointing out to chem- 
ical variations across the various velocity regimes at different 
positions. 
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Appendix A: Diagnostic diagrams 

The results of the excitation analysis, presented in Sect. |4] are 
here illustrated through the use of diagnostic plots, employing 
both line ratios and integrated intensities of the water lines. Each 
H2O emission component, analyzed in the paper, is shown sepa- 
rately: the two velocity components in the water emission from 
the R4 position, R4-LV and R4-HV, are shown in Figs. lA.ll and 
IA.2I respectively, while the water emission from the B2 position 
is presented in Fig. IA.3I 

In particular, for each emission component we report a set 
of nine plots to explore the parameter (Tkin, hht, A^onho-HiO, 
0) space and visualize the results given in Table |2] and 
Fig. |4] We considered the three most significant line ratios 
and one integrated line intensity. In particular, for the R4- 
LV and R4-HV components, the 557GHz/1097GHz (green), 
1670GHz/1097GHz (blue) and 1113GHz/1097GHz (red) ratios 
and the integrated intensity of the 1097 GHz line are reported 
in each panel. For the B2 component, the 557GHz/l 1 13GHz 
(green), 1670GHz/l 113GHz (blue) and 752GHz/l 1 13GHz (yel- 
low) ratios and the integrated intensity of the 1113 GHz line are 
shown. 

At R4, to separate the contribution from the two velocity 
components with different excitation conditions along the line 
of sight (R4-LV and R4-HV), we divided the emission range of 
the water lines: between and 20 km s"' for the low velocity 
component, R4-LV, and between 20 and 60 km s"' for the high 
velocity component, R4-HV. 

As an example. Fig. lA.ll for R4-LV shows that, for a given 
temperature {T) and emission size (0), at low H2O column den- 
sities the line ratios only depend on the H2 density, without any 
dependence on the H2O column density; while for higher H2O 
column densities there is a degeneracy between the H2O column 
density and the H2 density, i.e. the line ratios depend on the prod- 
uct of both quantities. However, different line ratios and the inte- 
grated intensity of one line allow to select only a limited region 
of the plane and therefore a limited range of physical conditions 
(i.e. in the case of R4-LV, all the line ratios and the integrated 
intensity of the selected line only intersect in the bottom-right 
panel of Fig. lA.lb . 

The following conclusions can be drawn from the inspection 
of the diagrams in Figs. lA.lllA.3l 




R4-LV 



N 




13 14 15 13 14 15 13 14 15 

log„N(o-H,0) (cm-^) log,„N(o-H,0) (cm-') log,„N(o-H,0) (cm-») 

Fig.A.l: Diagnostic diagram to explore the parameter (r^n, 
"Hi, A^oitho-HTO, 0) space that fits the observed data of the R4- 
LV component. The different panels are for three different val- 
ues of the emission size (0 = 5", 15" and 35") and three 
different temperatures (T - 150, 300 and 650 K). Three line 
ratios are reported in each panel: 557GHz/1097GHz (green), 
1670GHz/1097GHz (blue) and 1 1 13GHz/1097GHz (red); and 
the integrated intensity of the 1097 GHz line (cyan). An error of 
10% is assumed for each line intensity. The thick boxes in this 
diagram and in all the following diagrams mark the range of pa- 
rameters for which all the line ratios and line intensity intersect, 
corresponding to the physical conditions that are consistent with 
the observed data. 




13 14 15 13 14 15 13 14 15 

log,„N(o-HjO) (om-2) log„N(o-HjO) (om"') log,„N(o-HjO) (cm-^) 



Fig. A.2: Same as Fig.lAHfor R4-HV. 
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Fig. A.4: H2O line profile (black) decomposition in a triangular low-velocity component (red), derived by fitting the 752 GHz water 
line to the low velocity emission profile, and a high-velocity component (green), obtained by subtracting the low-velocity component 
from the total water line profile. 
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Fig. A. 3: Diagnostic diagram to explore the parameter (Tkin, n-a^, 
A^ortho-HiO, 0) space for the water emission at the B2 position. 
The different panels are for three different values of the emis- 
sion size (0 = 5", 20" and 35") and three different tempera- 
tures {T - 150, 450 and 650 K). Three line ratios are reported 
in each panel: 557GHz/1113GHz (green), 1670GHz/1113GHz 
(blue) and 752GHz/1113GHz (yellow); and the integrated inten- 
sity of the 1113 GHz line (cyan). An error of 10% is assumed 
for each line intensity. 



1 . The water emission at R4-LV is consistent with an extended 
component, with a very dense gas (larger than 10^ cm"^) and 
low column density (a few lO'"* cm"^). The temperature of 
this extended R4-LV component is high, about 600 K. 



2. R4-HV, at lower excitation than R4-LV, can be reproduced 
by a less extended emission (about 15-20 arcsec) with a 
range of temperatures from 150 K to 650 K and either a 
higher density (larger than 10^ cm"^) and lower column den- 
sity (less than lO''* cm"^) or a lower density (around lO*" 
cm"^) and higher column density (a few lO'"* cm"^), respec- 
tively. 

3. The emission at B2 appears to be consistent with an extended 
component (around 35 arcsec), with a wide range of temper- 
atures, high density (larger than 10^ cm"^) and low column 
density (around or smaller than 10'"^ cm"^). However, we 
know from the PACS map at 1670 GHz towards L1448 (see 
Fig-ID that the water emission comes from a more compact 
region. Therefore, the solution for an emission size about 20 
arcsec was also considered, although it does not reproduce 
one of the line ratios (752GHz/l 1 13GHz) for any of the dis- 
played temperatures. Thus to further constrain the physical 
parameters of the water emission at the B2 position, some as- 
sumptions need to be made. A s explained in Sect. |4l Spitzer 
mid-lR H2 observations from iGiannini et al.l (1201 ll) has be 
used in this context to constrain the temperature to ~ 450 K, 
assuming that water emission originates from the same gas 
emitting the mid-lR H2 lines. 

Finally, as mentioned in Sect. |4] we also explored how the 
derived results can be affected by the method used to separate 
the high and low velocity components in R4. For that, we tried to 
deconvolve the line profiles into two different components based 
on the line shape: R4-LV being identified with the triangular pro- 
file seen especially in the higher excitation lines; while R4-HV 
corresponding to the high velocity emission component that is 
overlaid on the triangular shape in the lower excitation water 
lines. Given the triangular line profile of the para-HaO 2\i - 2o2 
line at 752 GHz, we assumed that all the emission in this line 
is coming from R4-LV. We thus used this line as a template for 
separating the R4-LV contribution in each water line profile, as 
illustrated in Fig. IA.4I we scaled the 752 GHz line to fit the 
low velocity triangular contribution (red curve) in each water 
line profile (black curve) and then subtracted it from the rest of 
the line, to obtain the R4-HV component (green curve). In this 
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Fig. A. 5: Diagnostic diagram to explore the parameter (Tyn, rta^, 
A^ortho-H20, 0) space for the water emission at R4-LV, after di- 
viding the two velocity components as explained in the text and 
illustrated in Fig. IA.4I The different panels are for three differ- 
ent values of the emission size (0 - 5", 15" and 35") and three 
different temperatures (T = 150, 300 and 650 K). Three line 
ratios are reported in each panel: 557GHz/1097GHz (green), 
1670GHZ/1 097GHz (blue) and 1 1 13GHz/1097GHz (red); and 
the integrated intensity of the 1097 GHz line (cyan). An error of 
10% is assumed for each line intensity, to simplify the plots. 



way we decomposed every water line profile in the two contri- 
butions from the R4-LV and R4-HV gas components. For each 
water line, both derived spectra have thus been integrated over 
the whole emission range to derive the integrated intensities of 
the two components. 

We thus compared the derived line ratios and one integrated 
Hne intensity of R4-LV (Fig. lA.5l ) and R4-HV (Fig. lA.6l l with the 
models. The considered integrated line intensities for each com- 
ponent are the 1097 GHz line for the R4-LV and the 1670 GHz 
line for R4-HV. Moreover, for R4-HV, we show also the upper 
limits to the 752GHz/1670GHz line ratio (yellow), derived from 
the 3cr upper limit of the 752GHz integrated intensity, in the only 
three panels where it is significant and can be used to discrimi- 
nate between the different models. 

The plots clearly highlight the same results obtained in the 
previous section, meaning that the approach we use to separate 
the two velocity components at different excitation at the R4 po- 
sition does not affect the results of the analysis. This is mostly 
because the two velocity components are well separated in ve- 
locity (R4-LV peaks at about 10 km s"' and R4-HV at about 
25 km s"') and therefore each of them dominates the water emis- 
sion in the relative velocity range (0-20 km s"' for R4-LV and 
20-60 km s ' for R4-HV) we chose in Sect.|4] 



e= 5" B= 15" e= 35" 
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Fig. A. 6: Same as Fig. IA.5I for R4-HV. Two line ratios 
are reported in each panel: 557GHz/1670GHz (green) and 
1113GHz/1670GHz (red); and the integrated intensity of the 
1670 GHz line (cyan). Moreover the upper limits to the 
752GHz/1670GHz (yellow) line ratio is shown in the right pan- 
els, corresponding to = 35". 
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